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In the study of entanglement in a spin chain, people often consider the nearest-neighbor spins. 
The motivation is the prevailing role of the short range interactions in creating quantum correlation 
between the 1st neighbor (IN) spins. Here, we address the same question between farther neighbor 
spins. We consider the one-dimensional (ID) spin-1/2 XY model in a magnetic held. Using the 
fermionization approach, we diagonalize the Hamiltonian of the system. Then, we provide the 
analytical results for entanglement between the 2nd, 3rd and 4th neighbor (denoted as 2N, 3N, and 
4N respectively) spins. We find a magnetic entanglement that starts from a critical entangled-field 
(he) at zero temperature. The critical entangled-field depends on the distance between the spins. In 
addition to the analytical results, the mentioned phenomenon is confirmed by the numerical Lanczos 
calculations. By adding the temperature to the model, the magnetic entanglement remains stable 
up to a critical temperature, T c . Our results show that entanglement spreads step by step to farther 
neighbors in the spin chain by reducing temperature. At first, the IN spins are entangled and then 
further neighbors will be entangled respectively. T c depends on the value of the magnetic field and 
will be maximized at the quantum critical held. 

PACS numbers: 03.67.Bg; 03.67.Hk; 75.10.Pq 

I. INTRODUCTION 

The quantification and control of quantum entan¬ 
glement, for quantum information processing purposes 
through teleportation^”— and superdense coding^, have 
derived extensive experimental^”— and theoretical^ — re¬ 
searches in the recent decades. Physically interest¬ 
ing in this area, solid state systems also have been 
used as a topic of many studies for their potential ap¬ 
plications on transmitting a quantum state by using 
entanglement^ 2 . Particularly one-dimensional quantum 
spin systems comprise many nonclassical properties that 
make them interesting for studying spin entanglement. 

The magnetic behavior of these systems is explained 
through the isotropic Heisenberg model or its special 
cases such as the Ising, XY, and XXZ models. In 
addition to many compounds with common isotropic 
Heisenberg structure like La^GuO 4 ^, CuGeOj —, and 
LiCuVO/r — , there are also some experimental examples 
such as Fe(C 5 H 5 N0)e(Cl0 4 )2r^ and Gs^CoCl^ which 
are well described by Ising and XY models respectively. 

It has been found that the systems with XY interaction 
can be used as quantum dotsi£”—. Therefore, many stud¬ 
ies have been devoted to the exploration of important 
features of this modeli^r— . Also the fact that its eigen¬ 
values can be exactly solved through the Jordan-Wigner 
(JW) transformation makes this model very practical— 

The zero-temperature quantum behavior of the spin- 
1/2 isotropic antiferromagnetic XY chain shows that it is 
in the Luttinger Liquid (LL) phase and undergoes a quan¬ 
tum phase transition by applying an external magnetic 
field~. In the LL regime, the IN spins are entangled— 

Increasing the magnetic field reduces the entanglement 
measure until the spins become completely disentangled 


in the magnetic fields larger than the quantum critical 
field, h c . 

The thermal entanglement (TE) is of particular inter¬ 
est and demonstrates that the non-local correlations per¬ 
sist even in the thermodynamic limi t 26 ! 27 . Studying the 
TE, zero-temperature entanglement, and the relation be¬ 
tween the entanglement and the quantum phase transi¬ 
tions may lead to a relationship between the quantum 
information theory and the condensed matter physics. 

Since the TE can be experimentally verified via mea¬ 
surable macroscopic parameters^”—, many efforts have 
been dedicated to quantify it. Also, it has been recently 
demonstrated that an entangled pair of spins separated 
by several lattice spacing within an antiferromagnetic 
chain can be a suitable candidate for application in the 
quantum information processing— This work also high¬ 
lights the importance of studying the TE in different dis¬ 
tances. 

According to Ref. ;33j, TE of the IN spins decreases by 
increasing the temperature. There is a critical tempera¬ 
ture (T c ) where TE will be zero— It was found that the 
critical temperature is independent of the magnetic field. 
This property is not exclusive to the XY model, but it has 
been observed in other ID spin-1/2 systems as wel l 34 ! 35 . 
In addition to the IN spins, the zero-temperature entan¬ 
glement between two spins at arbitrary distances was also 
studied using analytical and numerical calculations in the 
spin-1/2 XY chain s 36 ! 37 . It has been shown that at zero 
temperature, the entanglement between the non-nearest 
neighbor spins becomes finite around the quantum crit¬ 
ical field. The entanglement remains finite even if the 
distance between the spin pair reaches the system size— 
The concurrence and the quantum discord for the 2N 
spins were also studied in the Ising and the XXZ spin 
chains— It was found that the quantum correlations 
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increases in the region close to the critical points. The 
scaling behavior of the quantum correlations between a 
spin pair at an farther distances in the XY chain has been 
investigated recentl y 39 ' 40 . It is shown that the quantum 
discord between the spin pairs positioned at distances 
further than two lattice spacing reveals a quantum phase 
transition. Moreover, the quantum discord may increase 
both with the temperature and the magnetic field in cer¬ 
tain regions of the parameter space. 

In this paper, we study an infinite ID spin-1/2 isotropic 
XY model in the presence of a magnetic field. Using the 
JW transformation, we find an analytical solution for 
the TE between the 2N, 3N, and 4N spin pairs in the 
thermodynamic limit. In absence of the magnetic field, 
the non-nearest neighbor pairs are not entangled. They 
remain unentangled until the magnetic field reaches the 
critical value hf which depends on the distance between 
the spins. We show that the hf is a function of h c . By 
further increasing the magnetic field, the entanglement 
between the 2N, 3N, and 4N spin pairs increases to a 
maximum value before settling to zero at the quantum 
critical field, h c . By taking the temperature into account, 
the concurrence within 2N, 3N, and 4N pairs decreases 
and reaches zero at a magnetic field dependent on crit¬ 
ical value T c . Note that this critical temperature was 
found to be field-independent for the IN spin a 33 ' 35 . T c 
is maximum at the quantum critical field value. For the 
magnetic fields larger than the quantum critical value, 
two critical temperatures can be found. 

This paper is organized as follows. In the next section 
we introduce the model and the Hamiltonian which is 
diagonalized by using the JW transformation. In section 
m we discuss our results on the thermal behavior of 
the concurrence between the 2N, 3N, and 4N spins. In 
section m a conclusion and the summary of the results 
is presented. 


II. THE MODEL 

The Hamiltonian of a ID spin-1/2 isotropic XY model 
in a external magnetic field is written as 

N N 

n = JX>iS'« + s”s‘V 1 )-)>£s-, (i) 

1=1 1=1 

where Sj is the spin-1/2 operator of the j-th site. J > 0 
denotes the antiferromagnetic exchange coupling and h is 
the magnetic field. We use the JW transformation to di¬ 
agonalize the Hamiltonian. It converts the spin operators 
into spinless fermion operators as below: 


s; 

= a t(e iw5: «W“‘), 

(2) 

s 7 

= ( e -^£.<W“! )aj, 

(3) 

Si 

t 1 

= a ) a o - 2 ' 

(4) 


By performing this transformation, the Hamiltonian is 
mapped onto the Hamiltonian of a ID non-interacting 
fermionic system, 

H f = ^y + J 1 + a A+ 1 ) ~ h Yl a ]aj.(5) 

3 3 

Using the Fourier transformation aj = - 7 = J2k e ~ lk ^ a ki 
the momentum space Hamiltonian is diagonalized as 

Hf = J2 £ (k)a{a k . ( 6 ) 

k 

s(k) is the dispersion relation: 

e(fc) = J cos(k) — h. (7) 

III. THERMAL ENTANGLEMENT 

We focus on the entanglement between two sites, which 
is quantified by the concurrence. As defined below the 
concurrence measures the non-local quantumness of the 
correlation s 6 ' 7 ' 41 . 


C(pij) = max {0, Ai - A 2 - A 3 - A 4 }, ( 8 ) 


where A^s are square roots of the eigenvalues of the prod¬ 
uct matrix R = pp where = (erf ® crJ)pT(crf ® erf) 
and p^ is the reduced density matrix related to any 
(i,j) pair in the chain. A^s are real and non-negative 
even though R is not necessarily Hermitian^. The re¬ 
duced density matrix contains all information about the 
involved spins. The density matrix in the standard basis, 
{| tt)> I t 4 -)> I 4 -t)j I 4 - 4 -)}> can be expressed as: 


pij — 


/ < p/p/ > < P/ S~ > < S rpj > < S"S- > \ 

< P/St > < P/P/ > < S"St > < S-p/ > 

< S+P/ > < S+S” > < p/p/ > < p/s- > 

\ < S+S/ > < s+p/ > < P/S+ > < p/p/ > 


The brackets symbolize the ground state and the ther¬ 
modynamic average values at zero and a finite temper¬ 
ature, respectively, pt = \ + S z ,pJ = \ — S z , and 
S+ = S 31 ± iS v . Following the symmetry properties of 
the Hamiltonian, the density matrix must be real and 
symmetrical^. It is found that only some elements of 
the density matrix are non-zer o 44 ' 45 : 


Pij — 


(K 0 0 0 \ 

0 Y+ Z* 0 0 

0 Zjj YP 0 

\ 0 0 0 XT.) 


( 10 ) 


By considering the spin pair (i,j) at the distance m from 
each other, j = i + m, the density matrix elements can 
be obtained as: 
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h 



obtained as follows: 

Zi,i+2 = /2-2/0/2+2A 2 , 

X+i+2 = fl ~ fl (13) 

Zi,i+3 = 4(/| ~ 2/o/l/2 + /|/l + /o/3 

— /1/3 + /1/2 — /0/3) + /3; 

= fl ~ fl (14) 

^i.i+4 = 8(/-f — 3/o/i /2 + 2/i/| + 2 /q/i/3 

+ flfl - fl - 2 / 0 / 1 / 2/3 + 2 hflh - Vlh 
+ flfl - / 0 / 2 / 3 2 - flh + 2 /oA 2 /4 - 2/ 2 / 2 / 4 
+ foflU) + 4(3/ 2 / 2 - 2 / 0 /I - 4 / 0 A /3 
+ 2AAA + 3/q/ 4 - 2/ 2 A + A/| - A 2 / 4 ) 

+ 2(2/iA — 3/oA + /|) + A) 

= fl - fl (15) 

where for the non-negative integer number n 

fn=7 tf e ~ ikn f^ dk - ( 16 ) 


FIG. 1: (color online). The concurrence between (a) 2N, (b) 
3N, and (c) 4N spins versus the applied magnetic field at zero 
temperature. The inset of (a) shows Lanczos results for the 
2N spins obtained from the finite chains with lengths IV = 
16, 20, 24 and exchange coupling J = 1. 



— < 77-2?T»2_|_ m >, 

Y + 

— < Tli( 1 77^_|_ m ) , 

y- 

< Hi) ? 


= < a\(l - 2a\ai)(l - 2a\ +1 a i+ i) 



V 

— <^1 Til 77/2_|_ m + , (11) 

where n* = ajaj is the occupation number operator. 
Therefore, the concurrence between this spin pair is: 


C(p) = max{0,2(|Z M+m | - X+ +rn X~ l+ J}. (12) 


The to = 1 case has been discussed in Ref. [13|. Gong 
and Su have derived an equation for a unique T c so that 
the IN spins are entangled at temperatures lower than 
T c . Here we study the concurrence of the 2N, 3N, and 
4N spins only. For these cases Z i}i+m and Xf i+m are 


/(*) = I+i is the Fermi distribution function, 
where /3 = and the Boltzmann constant Icb = 1. 
Eq. (fl3l) . Eq. (fbH) . and Eq. (fl5l) are obtained using the 
Wick’s theorem^. 

In a system with a fixed exchange coupling, these in¬ 
tegrals depend on two control parameters, i.e. tempera¬ 
ture (T) and magnetic field ( h ). A particular interesting 
case is study of the quantum correlation between the 2N, 
3N, and 4N spins at zero temperature. Fig[l] depicts the 
analytically calculated concurrence as a function of the 
magnetic field. Note that although the 2N, 3N, and 4N 
spin pairs follow the same trend in FigsQJa), (b), and (c) 
the concurrence measure decreases by increasing the dis¬ 
tance. Also none of the pairs is entangled in the absence 
of the magnetic field. They remain unentangled up to a 
critical entangled-field, hf. Fig|T] shows that hf for the 
4N spins is closer to the quantum critical field value than 
the other two spin pairs. It is notable that there is a re¬ 
lationship between hf and h c . Our results show that the 
critical entangled-field is 4/i c and approximately | h c and 
jijh c for the 2N, 3N, and 4N pairs respectively. Based on 
these results, one may propose a general relation for hf 
as a function of the distance and h c as hf' = , ,, h c . 

According to this relation and in agreement with the 
results from the finite size system^, it seems that all 
neighbors are entangled in the vicinity of the quantum 
critical field. Increasing the magnetic field above the hf 
induces entanglement between the non-nearest neighbor 
spins. This entanglement that is induced by the external 
magnetic field, is known as the ’’magnetic entanglement” 
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FIG. 2: (color online). The temperature dependence of the 
concurrence between (a) and (b) 2N, (c) and (d) 3N, and (e) 
and (f) 4N spins. The left-hand (right-hand) figures are for 
magnetic field values less (larger) than the quantum critical 
field. 


and was reported in finite size systems for the first time^. 
Here we explicitly showed that this phenomenon is also 
observed in the thermodynamic limit. Furthermore, the 
concurrence increases to a maximum value before falling 
to zero at the quantum critical field, h c . In the saturated 
ferromagnetic phase, h > h c , the non-nearest neighbor 
spins are not entangled similar to the IN spins. We also 
calculated the spin-spin correlation functions < > 

for a = x,y,z. We found that the spin-spin correla¬ 
tion functions along the x and y directions dominate the 
correlations along the z direction exactly at the critical 
entangled-field. We have also calculated the concurrence 
between the 2N spins in finite size chains by using the nu¬ 
merical Lanczos method. The numerical results for chain 
sizes N = 16, 20, 24 are plotted in the inset of Fig. CD (a). 
Although the effect of the level-crossing is seen in the fi¬ 
nite size results, there is a region where the 2N spins will 
be entangled by the magnetic field. This is in agreement 
with the exact analytical results. 

To obtain a better insight into the nature of the en¬ 
tanglement between the non-nearest neighbor spins, we 
consider the thermal behavior of the entanglement as a 
general case. Fig. [2] shows the temperature dependence 
of the concurrence between the 2N (Figs. [2] (a) and (b)), 
3N (Figs. [2] (c) and (d)), and 4N (Figs. [5] (e) and (f)) spins 


for different values of the magnetic field h> hf. Figs. [2] 
(a), (c), and (e) show the results for hf < h < (h c ). The 
2N, 3N, and 4N spins are entangled at zero temperature 
in this region. By increasing the temperature, the en¬ 
tanglement decreases and vanishes at a critical value T c . 
The existence of a critical temperature is also reported 
for concurrence between the IN spins in this model^. 
Here, the critical temperature is a function of the mag¬ 
netic field strength and increases when the magnetic field 
approaches the quantum critical value. It is interesting 
that the critical temperature (T c ) has a maximum exactly 
at h = h c = 1.0. On the other hand, the maximum value 
of the critical temperature for the 3N spins in this model 
is almost one-third of the maximum value of the critical 
temperature for the 2N spins and three times larger than 
the maximum value of the critical temperature for 4N 
spins. By comparing Figsj2](a), (c), and (e) we conclude 
that the entanglement is formed step by step from the 
2N, to the 3N and to farther neighbors in the spin chain 
by reducing the temperature. So the 4N spins are entan¬ 
gled at a lower temperature than the other two pairs. 

For the magnetic fields greater than the quantum crit¬ 
ical value (Figs. [2] (b), (d) and (f)) the 2N, 3N, and 4N 
spins are not entangled at zero temperature. By increas¬ 
ing the temperature from zero, the spins remain unentan¬ 
gled until the first critical temperature T Cl (h) is reached. 
Once the temperature exceeds T Cl , the concurrence re¬ 
gains and takes a maximum value, then decreases to zero 
at the second critical temperature T C2 (h). The existence 
of the second critical temperature is expected, since a suf¬ 
ficiently large thermal fluctuation will destroy all of the 
classical and quantum correlations. Both of the T Cl and 
T C2 values increase by increasing the external magnetic 
field. We have to mention that the same phenomenon has 
been observed for the concurrence between the IN spins 
in this region of the magnetic fields. The only difference 
is that the second critical temperature for the concur¬ 
rence between the IN spins is almost field-independent. 
Also note that the entangled region beyond the h c is very 
limited for the 3N and 4N spins and diminishes quickly 
after h = 1.1 in both cases. 

Fig. [3] depicts the T-h phase diagram for (a) 2N, (b) 
3N, and (c) 4N pairs. Also Fig. [3] (d) shows an overview 
of the entangled regions for the 2N, 3N, and 4N pairs. 
These diagrams signify that there is no entanglement at 
any temperature interval in the absence of a magnetic 
field. The borders of each entangled region determine 
critical temperatures in each magnetic field. It is obvious 
that T c i = 0 and T c 2 has an incremental behavior in 
hf < h < h c for all three pairs. For h > h c , however the 
increasing trend of T c i can be observed for the 2N, 3N, 
and 4N pairs. T c 2 shows a different behavior and T c 2 is 
almost independent of the magnetic field for the 2N spins, 
while it decreases with h for the 3N and 4N. According 
to Fig. [3] (d), one can easily conclude that increasing the 
magnetic field and decreasing the temperature leads to 
the induction of entanglement on the farther neighbor 
spins. Therefore, all pairs are entangled around the h c 
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FIG. 3: (color online). Schematic T-h phase diagram for (a) 
2N, (b) 3N, (c) 4N spin pairs. Panel (d) shows the overview 
of entangled region for all three pairs. 


at zero temperature. 

The complete phase diagram of the concurrence be¬ 
tween the 2N, 3N, and 4N spins in the spin-1/2 XY chains 
in a magnetic field is plotted in Figs. |4] (a), (b), and (c) 
using the exact analytical results. The concurrence re¬ 
duces by increasing the temperature in all of the men¬ 
tioned spin pairs. However, the concurrence magnitude 
and the thermal interval of the entangled region for the 
3N and 4N spins are considerably smaller than the 2N 
spins. 


found that the thermal entanglement becomes finite only 
after h crosses hf value. When h increases beyond the 
critical point h c , the entanglement is still non-zero and 
decays gradually which is different from the ground state. 
The entangled field region is smaller at higher tempera¬ 
tures. Furthermore, unlike the IN spins entanglement, 
there is not any unique critical temperature for the non- 
nearest neighbor spins. 
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IV. CONCLUSION 

In this work we proposed a generalizable method for 
calculating the TE between any arbitrary pair of spins in 
low-dimensional spin chains via the analytical JW trans¬ 
formation. We exactly resolved the concurrence between 
the 2N, 3N, and 4N spins in the spin-1/2 XY chains in 
the presence of a magnetic field. We found that all of 
the spin pairs are entangled at zero temperature where 
hf < h < h c . We also showed that the critical entangled- 
field, hf, is a fraction of the quantum critical value. The 
critical entangled-field is h c /2 for 2N spins and ~ 4h c /5 
and 9/i c /10 for the 3N and 4N spins respectively. Fur¬ 
thermore, we studied the temperature dependence of the 
concurrence between the 2N, 3N, and 4N spins. It is 



FIG. 4: (color online). The curved surface of the concurrence 
between (a) 2N, (b) 3N and (c) 4N spins as a function of 
the temperature and the applied magnetic field for exchange 
coupling J — 1. 
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